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^~~^ Abstract: Formalism to calculate the hydrodynamic fluctuations by applying the On- 

r~i sager theory to the relativistic Navier-Stokes equation is already known. In this work, 

Oj we calculate hydrodynamic-fluctuations in the framework of the causal hydrodynamics of 

^ Miiller, Israel and Stewart and the other related approaches. We show that expressions 

O for the Onsager-coefficients and the correlation-functions have forms similar to the ones 

obtained by using the relativistic Navier-Stokes equation. However, spatio-temporal evo- 

^ lution of the correlation functions obtained using MIS and the other causal theories can be 

J> significantly different from the correlation functions obtained using the Navier-Stokes equa- 

r^ tion. Finally, as an illustrative example, we numerically evaluate the correlation-functions 

00 using the one dimensional expanding boost-invariant (Bjorken) flows and compare the 

correlation-functions obtained using the various hydrodynamic approaches. 
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1. Introduction 

A study of fluctuations in continuous media is of great interest in physics and it can 
provide a link between the macroscopic and microscopic points of view. A macroscopic 
theory such as hydrodynamics provides a simplest possible description of a complicated 
many-body system in terms of space-time evolution of the mean or averaged quantities like 
energy density, pressure, averaged flow velocity etc. The fluctuation theory studies small 
deviations from the mean behavior and it can help in calculating correlation functions for 
the macroscopic variables [1,2]. In context of relativistic hydrodynamics, results of the 
fluctuation-dissipation theorem has been studied in Refs. [3,4]. In Ref. [3] authors have 
studied the fluctuation in the contexts of general-relativistic Navier-Stokes theory. A more 
general framework of hydrodynamics described as the divergence type theory(DTT) [5] was 
considered in Ref. [4]. It ought to be noted that recently in an interesting work in Ref. [6], 
the authors have applied the results of fluctuation-dissipation theorem to the relativistic 
Navier-Stokes theory of hydrodynamics and calculated the two particles correlators for 



the one-dimensional hydrodynamics (Bjorken) flow relevant for the relativistic heavy-ion 
collision experiments at RHIC and LHC. The authors obtained several analytical results 
for two particle correlation functions. Further, in Ref. [7], the authors have studied the 
effect of thermal conductivity on the correlation function using the Bjorken-flow. It should 
be noted here that it is well-known that relativistic Navier-Stokes theory exhibit acausal 
behavior and it can give rise to unphysical instabilities [8] . However the causality can be 
restored if the terms with higher orders are included in the hydrodynamics as indicated 
by the Maxwell-Cattaneo law [9]. Indeed these issues do not arise in the second-order 
causal hydrodynamics theory developed by Miiller, Israel and Stewart (MIS) [10]. The 
Navier-Stokes equations can be derived from the second law of thermodynamics d^^S^ > 0, 
where S^ denotes the equilibrium entropy current [2]. However, it is not necessary for an 
out-of-equilibrium fluid to have an equilibrium entropy current [9]. In MIS hydrodynamics 
out-of-equilibrium current can have contributions from dissipative processes like the effect 
of viscosity and the heat conduction. This has an interesting analogy with the irreversible 
thermodynamics [11, 12]. Further the MIS hydrodynamics has been extensively applied 
to study the relativistic heavy-ion collisions [9, 13-15] and also in cosmology [16]. Later 
this formalism was extended to include the effect of third order terms in the gradient 
expansion [17]. Recently, it has been shown that the derivation of the MIS equations from 
the underlying kinetic equation may not be unique, there may exist a more general set 
of hydrodynamic equations which may allow one to obtain MIS equations as a special 
case [18,19]. 

Finally, it should be mentioned here that although the divergent type theory (DTT) 
of relativistic fluid of Geroch-Lindblom [5] allows for a consistent proof of causality and 
stability of its solutions, it is far from direct thermodynamic intuition. Moreover, the 
connection between the DTTs and MIS or other causal hydrodynamics theories is not yet 
clearly established. 

In this work we apply the fluctuation-dissipation theorem to MIS equations and also to 
the hydrodynamics models developed by Denicol, Koide and Rischke (DKR) [19], Jaiswal, 
Bhalerao and Pal (JBP) [18] and other models based on MIS approach [10,17,20] Further, 
we apply the results obtained to study the hydrodynamical evolution using 1-1-0 dimensional 
Bjorken flow. In particular, we calculate the correlators using the Onsager coefficients for 
various relativistic hydrodynamical theories. We believe that such study may be important 
in choosing the appropriate version of the causal hydrodynamics from the variety of the 
options available. 

2. Fluctuations and correlations in Hydrodynamics 

In thermodynamic equilibrium entropy of the system S which is function of the additive 
quantities Xj. becomes maximum. In equilibrium, S satisfies the condition Xj. = — ^- = 0. 
However, when the system is slightly away from the equilibrium the generalized forces 
Xf. / and -^ = —^^j^Xj^ + ^j, the summation convention is implied, describes the flux 
associated with the quantity Xj, here ^j are the random forces or the noise term and ^^j^ are 
the Onsager coefficients. The Onsager reciprocity relations imply that 7j^ = 7^^. In this 



phenomenological theory time rate of change of the total entropy, -^ can be given by, 

^ = -^X„ (2.1) 

dt dt ^ ^ 

which can also be written as, 

-17 = lik^kXi - (iXi. (2.2) 

Noises are correlated and the correlation is given by the formula, 

{Uh)Ut2)) = {l^k + IkMh - h)- (2.3) 

The correlation functions can be found once 7j,- are known [2,6]. In order to find out 
7jj, one needs to know the -^ for the hydro dynamical theory and identify the generalized 
forces and fluxes. The expression for entropy rate -^ can be found either by using equations 
of hydrodynamics together with the thermodynamic relations or from kinetic theory [10]. 
In our work we are using the former approach. 

We start with the expressions for the energy momentum tensor T^'^ and the current- 
density J^ for a viscous fluid, 

j^ = ubu^ + y^ + I", (2.5) 

where, Tj^ = eu^u^ — pA^'^ is the ideal part of the energy momentum tensor with 
A'^'^ = g^'^ — u^u'^ , and e, p, u^ are the local energy density, pressure and fluid flow 
four-velocity respectively. It is to be noted that u^u^ = 1 and g^^j = diag{+,—^—,—). 
AT^*^ = AT,^.^ + AT^^;„, with ATZ = ^^" - A^"n and AT^^^^, = W^u"" + Wu^, is the 
dissipative part of the energy momentum tensor and S^^ is the stochastic term arising 
due to the local thermal fluctuations [6]. Similarly, v^ and l'^ denote the dissipative(non- 
equilibrium) and the stochastic terms in baryon current density respectively, n^ denotes 
the net density of baryon number in local rest frame. W'^ = q^ + hv^ is the energy flow in 
local rest frame, h = {e + p)/n is the enthalpy per particle and u^ = A^^ J^u is the baryon 
number flow in the local rest frame. For the dissipative fluxes one can always require the 
following relations to hold u^tt^^ = 0, tt^ = 0, A^j^tt'^'^ = 0, u^W^ = 0, u^u^ = 0, u^q^ = 0. 

The relevant conservation equations for the hydrodynamics can be written as, 

a^ J^ = DriB + ubV^u^ + d^^v^ = 0, (2.6) 

u^d^T^" = De + {e + p + Ii)V ^u^ - tt^^V^^'u'^ + V ^,W^ - 2W^Du^ = 0, (2.7) 
A'^d^T^" = {e + p + n)Dn" - V"(p + H) + A^'^V^vr,, - ^"'^Dn, 

+A'"'DW^ + IW^'^Vy^ = 0, (2.8) 



where, D = u'^d^ and V^ = A^^ d^. There are only 5 equations written above and 14 
unknowns n^, e, 11, VF^, t:^^ and u^ . Therefore, 9 additional equations for dissipative 
fluxes are required to obtain the hydro dynamical solution. 



There are two popular choices for u'^: In Landau-Lifshitz frame u^ is parahel to the 
energy-flow and W'^ = which implies that q'^ = —hv^. Another choice is the Eckart-frame, 
where u'^ can be parallel to J^/i and v^ = {) and this would imply W'^ = q^ . Now we shall 
obtain the 9 additional equations and fluctuation correlations in Landau-Lifshitz and Eckart 
frame respectively. 

2.1 Equations for dissipative fluxes in Landau-Lifshitz frame and fluctuation 
correlations in MIS 

In order to derive the 9 additional equations one needs the expression for the out-of- 
equilibrium entropy four-current. In Landau-Lifshitz frame the expression for the non- 
equihbrium entropy four-current is given in literature [10, 14, 15] and is as follows. 



5^ = su^ - ^u^ 



« rr2 « i^ I fl u\\u'' a^Iiq^ a^-n^^q^ 

PqH - /3iq^q + h-KyxT^ j — ^^ + . (2.9) 



The divergence of the non-equilibrium entropy four-current (Eq.(2.9)) using Eqs.(2.6-2.7) 
and the thermodynamic relations de = Tds — fidn and e + p = Ts — /xn can be written as 
follows, 



Td^S'' 



n 



d^u^ + /3oII + ^Td^ (^u^ ) n + aoV^q^ 



-h-'TV, f 



PiQ,, 



1 



Td,, 



T 



% - aiSj.TrJ^ + oif^d^Ii 



+ VT' 



fiv 



(^f,u - P2'^,,u - ^Tdx [ -^u^ ) TTf,^ + aiV(^g^) 



1. 



T 



(2.10) 



where we have used the notation F = DF. According to the second law of thermodynamics 
we must have Td^S^ > 0. This inequality will be satisfied if 11, q^^ and -k^^ satisfy the 
following equations. 



-h'^TV, 



^ 






1 



(^ijiv - ^T^fiu - ^Tdx [ Y^^ 



/?: 



T^fiu + aiV(^g^) 



n 

"C' 

XT' 
2r] ' 



(2.11) 
(2.12) 
(2.13) 



Here we note that the term with factor 1/2 on the left hand side of equations (2.11- 
2.13) are generally ignored by arguing that the gradient of thermodynamic quantities are 
small [15,21]. Thus, Eqs.(2.11-2.13) can be written as. 



-CV„n" - /ngV^g^, 



rnH + n 

Tqqf, + qi^ = -XT^h'^V^, [tJ ^ ^gnV^^n - /g^V^vr)^, 



(2.14) 
(2.15) 
(2.16) 



Henceforth, we call Eqs.(2.14-2.16) as MIS equations. Here, Tji = (/]q, t^ = \Tj3i, r^ = 2?7/32 
are identified as the relaxation times and /ng = Cc^O; ^gn = ATag, /q,r = \Ta-i l^g = 2r]ai 
as coupling constants. These 9 equations for the dissipative fluxes combined with the 
Eqs. (2.6-2.8) form complete set of the hydrodynamic equations. Note that the limit 
Tu^Tqj'^m ^-UqjlqUj ^qiri Kq " ^ is the first Order limit which correspond to the Navier-Stokes 
case. 

Eq.(2.10) can be written as. 



Td^S^ 



c 



XT 



+ 



III/ 



III/ 



2r] 



> 0. 



(2.17) 



Here, q^q^ < [13]. Now using the identity A'^'^A^^, = 3 and the condition A^^tt^'^ = 0, 
one can write Eq.(2.17) as. 



^ T \2r] 



A^.n 



3C 



XT^' 



Upon integrating over the whole volume Eq.(2.18) can be written as. 



dS 



^3 
a X 



AT'"' /vr 
2r? 



3C 






(2.18) 



(2.19) 



Following identification between the phenomenological variables (i^, X2) and the hydrody- 
namical variables can be made [6], 

X, ^ ATZ , X2 ^ g^ (2.20) 

A comparison of Eq.(2.19) with the phenomenological equation Eq.(2.1) will give, 

1 /vr„ 



X2 



'III/ 
T \2r] 

"^".AV. 



'"' > AF, 



3C 



(2.21) 



XT' 



Now neglecting the stochastic term in Eq.(2.2) and comparing it with Eq.(2.19) one 
can get. 



X, = -ATZ, 

I22X2 = -q^, 

7i2 = 721 = 



(2.22) 
(2.23) 
(2.24) 



The coefficients 7^2 and 721 are zero, because dissipative fluxes due to heat and viscosity 
are treated to be independent. Coefficients ^n and 722 are the rank- four tensor and they 
can be parameterized as follows. 



7ii 



j^^liuap ^ ^^liu^ap- 



1 

Ay' 



CA 



111/ 



722 



Ay 



(2.25) 



where, A 



/ii/alB 



/la ^^ufS 1 



A''"A 



A'"^A" . Now using Eqs.(2.22, 2.23) one can determine the 



coefficients A = 2r]T, B = QT and C = —XT . Thus one can write, 



711 = 2T 



fia ..ufS 



r/A'^^A 



o J Zi 



1 

Ay' 



722 



AT^A'''' 
AF ■ 



(2.26) 



Prom above expression of 711 one can see that there is a additive contribution of shear and 
bulk viscosity i.e one can write it as ^x\ = (7ii)?7 + (7ii)c- 

Now following Eq.(2.2), the correlation functions can be written as. 



(>5r.(^i)5.t(^2)) = 2r 



(/''(xi)r(x2)) 

(5r.(rri)r(x2)) 



r/(A^"A^^ + A'^'^A''") + (C - -?7)A^'^A°^ 

O 

n2 \^V 



2AT^A'"^(5(2;i-X2), 



0. 



5(xi-X2), (2.27) 

(2.28) 
(2.29) 



These are the stochastic or noise auto-correlation functions for the MIS hydrodynamics in 
the Landau-Lifshitz frame. 

2.2 Equations for dissipative fluxes in Eckart frame and fluctuation correlations 
in MIS 

In Eckart frame the expression for the entropy four-current [10, 13, 14] can be written as. 



,M 



5^ = s-" + ^ - (/5on - PiQ^q +/32vr^A7r j^ J^ ^ f ' 



(2.30) 



Note that here coefficients /3o , (^2 are same as in Landau-Lifshitz case while the coefficients 
di and /3i are given as, Oj = aj -|- ^^ and Pi = I5i + jj|— , where aj. Pi are the coefficients 
in Landau-Lifshitz frame. Next, the divergence of the non-equilibrium entropy four-current 
(Eq.(2.30)) using Eqs.(2.6-2.7) and the thermodynamic relations de = Tds — fxdn and 
€ + p = Ts — fin can be written as follows. 



Td^S^ 



n 



d^u^ + /3ori + l-Td^ (^u^ ) n + aoV^g'' 



V^lnT -iif,- Piq^ 



Td„ -^u q^ - aid^Tr'i^ + aod^U 



+ TT' 



/lU 



^iiu - h-^fMu - 7.Tdx ( -^u'' ) vr^j, + aiV(^g^) 



In order to have Td^S^ > we must have. 



(2.31) 






d.u^ + /3on + -Td. [^unn + a^V.q^ 



V^lnT -Uf,- Piq^ 



a 



flU 



-Td^ i Y^ J q^ - aid^n"^ + Ood^U 



n 

"C' 

AT' 
2r] ' 



(2.32) 
(2.33) 
(2.34) 



Thus one can write Eq.(2.31) as, 



Td„S^' 



tt2 a* a*!' 



c 



AT 



+ 



liv 



2r] 



>0, 



which can also be written as, 



(2.35) 






1^ 3C 2XT' 

Upon integrating over whole volume Eq.(2.36) can be written as. 



vr,, 



A^.n 



1 



Xuuqn+u q^) 



dS 
dt 






'^IJ.U ^t,u^ 



1 



2r/ 



^-l^[u,% + u,q. 



(2.36) 



(2.37) 



which can be rearranged in the following form. 



dS_ 
Itt 



d X 



^Z f^,u A^,n 



T 



27? 3C 



+ 






T^fjLv A^^n 



^ - -^^) - ^ ^^'^% + ^.^^) 



In this case also one can make the identifications as before, 

±1 -?■ AT^-^ , ±2 -^ ^Tf^eaf 

The comparison between Eqs. (2.38) and (2.1) will give, 

1 /7r„ 



X, 



X. 



' flU 



T \2ri 



1 
T 



IT 



fiu 



Au^n\ 
A^.n 



(2.38) 
(2.39) 

(2.40) 



2rj 3C 



1 



2 XT 



Uyqn + U^qu 



AV. 



Again neglecting the stochastic term in Eq.(2.2) and comparing it with Eq.(2.38) one can 
get. 



711^1 = 

722-^^2 = 

7i2 = 721 = 0. 



\rpfJ,y 



Now parameterizing ^n and 722 as follows. 



7ii 



AA 



/.luafS 



+ BA'^'A 



liv \alz 



1 

Ay 



, 722 



iA^^n^n'^ + BA^^n^'n" 



AV 



(2.41) 
(2.42) 
(2.43) 



(2.44) 



since we know the form of (X|,X2) and (AT^^^, AT^g'^J therefore using the above 
parameterizations of ^n and 722 and Eqs. (2. 41-2. 42), one can determine the coefficients 
A = 2r]T, B = QT and A = B = —2XT . Therefore the Onsager coefficients are. 



711 = 2T 



(riA^'^A"'^ - ^A'^^A'*'^') + ^CA'^'^A'^^ 



1 



722 



-2XT^ 



.A'a„,i',,/3 



Alia up, AUpil: 
^ U U + A ^ WU 



1 

Ay' 



(2.45) 
(2.46) 



Thus one can write the correlation functions using Eq.(2.3) as, 






2T 



^^^m^yp + A^^A''") + (C - -r?)A^^A"'^ 



-2Ar' 



+ A U U 



A^ u u + A u u + A'^ u u 
S{xi -X2), 



{sZi^i)s^Li^2)) = 0. 



5{x^-X2), (2.47) 

(2.48) 
(2.49) 



The form of these correlations is very similar to the correlations obtained for the 
relativistic Navier-Stokes case [6] . However, the evolution of the correlations can be very 
different as demonstrated later. 

2.3 Equations for dissipative fluxes in Landau-Lifshitz frame and fluctuation 
correlations for other Hydrodynamic models 



In this section we consider some of the interesting alternate approaches to the causal MIS 
hydrodynamics and some of its extensions. 

2.3.1 Third order hydrodynamics 

In the Ref. [17] third order corrections to the MIS hydrodynamics was considered when 
the effect of bulk- viscosity and heat-flux were absent. In this case expression for the non- 
equilibrium entropy four-current can be written as. 



P2 TjaB a , P2 a Bit u 



5^ = su'^ — ^TTnB^^'^u'^ + a^irvr^flVr^Tr'^'^u^, (2.50) 

where, q is a new dimensionless coefficient and it is assumed to be a constant. The last term 
on the right hand side of the above equation represents the third order correction to the 
equation of entropy. In order to fulfill the requirement of maximal entropy at equilibrium, 
the third order term must satisfy the condition a^ir^a'K'^'K u^ < 0. Divergence of the 
entropy four-current can be written as. 



d„S^ 



1 



TTajiCF 



ad,. 



al3 



VT, 



al3 



an n 



/3: 



2 n 



2T 



P2 ii\ a I3a , ^ a P2 a ■ /3a ^ r. 



1^2 .ajB 

T 



(2.51) 



Here, the Knudsen number(=r^^) is required to satisfy the condition t^9 <^ 1 for the 

validity of hydrodynamic approach. For the condition Td^S^ > to be satisfied one must 

have, 

1 

(2.52) 



'^>'' 2rjT 



III/ 

TT vr 



which implies that the form of shear viscous tensor vr should be given by 



fj,i/j 

alB 



vr 
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2'qT 
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;0" 
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2 . a/3 
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+ ad,, 
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(2.53) 



Now starting from Eq.(2.52) and following similar prescription to determine the Onsager 
coefficient as done in second-order MIS hydrodynamics one can write 

1 



X = TT' 



111/ 



X 



2t]T 



TT^^AV, 



and the Onsager coefficient, 

j = 2r]T 
The viscous correlation function can be written as 
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AF 



{S^-s{^i)St{x,)) = 2T 



r/(A^°A^^ + A'^^A 



fil3 ^^ua^ 



'-T^A^'A 



al3 



(2.54) 



(2.55) 



5{xi-X2). (2.56) 



One can notice that this expression is same as the one obtained using the second-order 
theory with 11 = 0. 

2.3.2 JBP hydrodynamics 

In the Ref [18] the authors have constructed the expression for the entropy four-current S^ 
generalized from the Boltzmann's H-function and find out the expression for its divergence 
as, 

n r. 



d,S^ 
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-(^annllV^Q + xa-^nT^lyuOL + Xann'^'iiK 
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fiu 
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flV 
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(2.57) 



where 6 = d^u^. The second law of thermodynamics Td^S^ > is guaranteed to be 
satisfied if we have, 
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therefore, vr, n and vr^ should satisfy the following equations. 






0-/.1/ - hT^t^u - /^Tvn^T^fiu - aiV(/,nj,) - xa,rn^(/.V^)a - xannn{f,u,y) 
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, (2.59) 
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where A, C, ^ > are the transport coefficients called coefficient of charge conductivity, bulk 
viscosity and shear viscosity respectively. Coefficients Oi, Pi^axY^ Pxx ^-re the additional 



transport coefficients and tlie parameters '0; X witli ip = 1 — ip and x = ^ — X included tlie 
contributions due to tlie cross terms of 11 and ii^'^ with n^ . 

Onsager coefficients in this case too, can be obtained using the parameterization[see 

Eq.(2.25)], 

1 



711 = 2r 



^A'^«A"'' + 1(^_^^)A^^A"^ 



Ay' 



(2.62) 



722 = ^^- (2.63) 



It should be noted that in Ref. [18] the authors have used ^ in Eq.(2.57) instead of ^ 
and therefore Onsager coefficient differ by factor T (see for example, Eqs.(2.26 and 2.63)). 
The correlation functions can be written as. 



(5r.(^i)c^(^2)) = 2T 



^^^m^^P + A^^A''") + (C - -r/)A^''A"'' 



5{x^-x^), (2.64) 



(/^(xi)r(x2)) = -2\T^^''5{x^ - X2), (2.65) 

(5„7,(xi)r(x2)) = 0. (2.66) 

2.3.3 DKR hydrodynamics 

In Ref. [19], it was demonstrated that derivation of relativistic viscous hydrodynamic 
equation from the 14-moment method done by Israel and Stewart may not be unique. In 
Ref. [19] , authors obtained relativistic dissipative hydrodynamic equations for the dissipative 
fluxes as, 

ti = ---M- ^liu^O + An.vr^V^,, (2.67) 



^IM-; = -— + 2p^a^'' + 27r^'^a;"^" - S^^t^^'O - r^^^'J^a'^'^ + A^nna^^ (2.68) 



M"") - JL I 9« ^/^'^ _L 9^<^,,'')" _ R ^''"a _ ^ ^(i^^")'^ _L \ „u^^^'' 



It should be noted that Eq.(2.68) contains vorticity term cj. This arises due to kinetic 
theory considerations and such term does not appear in the MIS hydrodynamics [13]. Note 
that in writing the above equations we have considered the fluid with no net baryon number. 
Thus the Eq.(2.7) with no net baryon number can be written as, 

9„(sun = :i:^ - 5^. (2.69) 

^.W^-l-^-'-m^ + ^I^^. (2.70) 

Pu Pum Pu Pn 

'"'^^ " 2/3, ^ 2/3,r, /3, ^ 2/3, ^ 2/3, 2/3, ' ^ ' ^ 

Assuming pressure gradient to be small and neglecting the higher (third) order term it is 
easy to show that -k^^^' = vr'^^. 

Neglecting the bulk viscosity and using Eq.(2.71) in Eq.(2.69) one gets. 
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where, S^ is the non-equihbrium entropy current for DKR hydrodynamics and has the 
form, 



S^ 



su 



aS a 



a/3 



'a^ 



a^ /3 



vr TTnRU T^^-n -EnRU 



2P^T 



TTTT" " Q/: 



4/3,r 



a/3 u a/^ H 



u a A 

4/3.T 

/i a A \ 



6/3,r 



2/3,T 2/3,T 



(2.73) 



Since Tr^*^ << J'^'^^ so 0{{tt^'^) ) terms will be smaller and its divergence will be much 
smaller. So on right hand side of Eq.(2.72) one can drop the divergence terms. Similarly, 
one can neglect the terms having divergence of four- velocity multiplied by 0{{'k^'^) ). Thus 
we can write Eq.(2.72) as. 



dS f 3 vr"'^ 



1 



2/3,r, 



(2.74) 



A comparison of the above expression with the phenomenological equation (Eq.(2.1)) yields. 



0/3 TA 



1 
f 



1 



Zp^T^ 



^aP^V. 



(2.75) 



Again by comparing Eq.(2.75) with Eq.(2.2)(when ^ = 0) one can get. 



(2.76) 



By considering 7 to be of the form. 



7 



^^/.^a/3 ^ ^^/.^^a/3- 



Ay' 



(2.77) 



and using Eq.(2.77) one finds that A = y — — — s- and i? = 0. Thus the viscous correlations 



2/3„ 



are. 



{SZi^,)St{x^)) = 2mr. (^A^'^A^^ + A^^A'^'^ - ^A'^^A"^) 5{x, - x,). (2.78) 



11 



2.3.4 Conformal viscous hydrodynamics 

The entropy current for the conformal hydrodynamics [20] can be written as, 



S^ 



su 



4riT 



n.«n"^^^ 



(2.79) 



One can easily find the following expression for the Onsager coefficient and the correlation 
function, 



7 = 27?T 



^Iia^ul3 



1 



liu xaP 



lA^'A 



{Svis{xi)S^is{^2)) = "^TlT 



^P I A^^l^ A'^°'\ 



^^t^a^uu ^ ^^,p^ 



1 



a/3 



(2.80) 
S{xi-X2). (2.81) 



3. Calculation of correlation functions in boost-invariant Hydrodynamics 

As an example we apply the results obtained in the previous sections to the relativistic 
heavy-ion collisions for the Bjorken scaling solution. 

According to Bjorken scenario in heavy ion collisions, the reaction volume is strongly 
expanded in the longitudinal direction, i.e along the collision axis(z-axis). So one can 
assume that there is no transverse flow. Thus one can describe flow in 1 + dimension [22]. 
It is useful to introduce the light cone variable y and proper time r which are defined by. 



T = \l t — z and y = arctanh(z/t) = -ln( ). 

The partial derivatives in time and space can be expressed as. 



cosh y — sinh y 
— sinh y cosh y 



Or 

y y. 



(3.1) 



(3.2) 



The flow velocity under the scaling assumption can be written as, u^ = 7(1, 0,0,^^) = 
(-,0,0,-) = (coshy, 0,0, sinh y). We consider only longitudinal flow fluctuations and 
parameterize the flow velocity [23] as. 



(cosh 6*, sinh/ 



(3.3) 



where 6 = y + 56{y, r) and 66(y, r) are the fluctuations in the longitudinal flow. In Bjorken 
scaling limit, 9 = y. With this parameterization and using the transformation of derivatives 
one can introduce the operators D, V such that. 



cosh(^ — y) sinh(^ — y) 
sinh(0 — y) cosh(0 — y) 



dr 

y y. 



(3.4) 



JL 

^M dr 

Since S^^ satisfies the condition. 



In the scaling limit, D = u^d^ 



dj- and 3„u^ = V^ 



u^S^'' = 0. 



One can write S^'^ as [6], 



S^'' = wiT)fiy,T)A^'', 



(3.5) 
(3.6) 
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where, w = e -\- p = Ts and / is a dimensionless quantity which satisfy (/) = 0, where 
() denotes the 'average value'. In heavy-ion collision experiments at LHC or RHIC a 
baryon free quark-gluon plasma is expected to be produced, therefore q^ = 0. Thus only 
viscous-correlations are of interest, which for MIS, JBP, third order and conformal case can 
be written as. 



(/(yi,n)/(y2,T-2)) = - — J^ 

AtiW (ji) 



3^(n) + C(n) 



^(n -^2)^(^1-2/2), (3.7) 



where 5{xi — X2)Transverse is replaced by effective transverse area A of colliding nuclei. Note 
that these correlations are same as that obtained by authors in Ref. [6] for Navier-Stokes 
case. Similarly, for DKR case one can write the viscous correlations as, 

(/(yi,n)/(y2,r2)) = —-^j-—^SiT,-T2)diyi-y2). (3.8) 

Atiw (ti) 

Defining 77 = /3^r^(see Ref. [19]) and neglecting the bulk viscosity for the correlation 
functions, one can rewrite the correlations for all the models of hydrodynamics as, 

(/(yi,n)/(y2,T2)) = ^^^<^(ri - T2)S{y, - y^), (3.9) 

where, 

^ ' 3riu;(ri) \s{ti) J ^j^^ 

Here, subscript [E] denotes the particular type of hydrodynamics model considered from 

the set of hydrodynamics models, for example [E] = [MIS, JBP, DKR, NS ]. 

It is useful to study the correlation function normalized by the initial value of the correlation 

obtained using the Navier-Stokes theory i.e. C(t) = — j- — ^ -^ — where, Tq is the 

•^ {to)ns^{to)ns 
initial-time for the hydrodynamics. C{t) can be written as, 

(iq\ (vil)\ "'(^) 

., \-lJ y^Jr) J [E] ^(r,) 

(^(V[E] = 7^(^yx • (3.11) 

Further, we neglect the effect of bulk-viscosity by considering the initial temperature Tj to 
be much larger than the critical temperature, T^ = 0.190 GeV. Now, in the Landau-Lifshitz 
frame, the energy and the momentum conservation laws are given by, 

u^d^T^" = De + {e + p)Ve - tt^^V^^u''^ - S^^¥^u^ = 0, (3.12) 

^^^d^T^"" = (e + p)Z)u" - V> + A^'^V'vr^^ - vr^'^Dn^ + A^'^a"^^^ = 0, (3.13) 

where, tt is the shear stress tensor and will be different for different models of hydrody- 
namics(e.g, MIS,JBP,DKR,NS etc.). 

In the scaling limit 9 = y, D = u^d^ = d^, d^u^ = V^ = - and shear stress tensor -k^^ 
has only diagonal components i.e ti^'^ = diag{0,Tr/2,TT/2, — vr). 
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Therefore, one can write the above equations as [19], 

9.e = -(^±^ + ^. (3.14) 

r T 

Here, the noise term is considered to be smaher than the background quantities. 

Equation for vr in the scahng hmit, for MIS, DKR and JBP hydrodynamics can be 
written as, 

9,7r + - = /3,;i-A-. (3.15) 

First we consider MIS case. In this hterature on relativistic heavy-ion colhsions A = case 
is being considered [15,21]. The coefficients /?^ and r^ can be written as. 

For JBP models parameters /3^, r^^ and A are as follows, 

2p _i 5ap 
/3,r = y, T,r =-Y'^ = 4/3> (3-17) 

where, a is the total cross-section [19] and it is assumed to be independent of energy [17,24]. 

In the massless gas limit the expressions for the coefficients f3^ and t^ are same in both 

JBP and MIS dynamics [18]. For DKR hydrodynamics, the parameters /3^, r^ [19] and 

A [25] are, 

4p _i 3crp 1 38 

Ptt = y > '^TT ^ 5 7^' 3^'^'^+ 'T'^ ^ 21' (3.18) 

In what follows we consider the ideal equation of state, e = 3p and the pressure is given by 

2 4 

the bag model, p = |qT . Further, we consider the initial values for temperature Tj = 0.310 
GeV and viscous stress vr = for the causal hydrodynamics and numerically solve Eqs. 
(3.14,3.15) for evaluating the correlations (3.11) in case of MIS, JBP, DKR hydrodynamics. 
However, for the Navier-Stokes hydrodynamics one needs to solve only Eq.(3.14) with same 
value of initial temperature and the viscous stress is given by, 

7r = r?^. (3.19) 

The results of the numerical work are presented in the following section. 

4. Results and Discussions 

We have studied fluctuations in various models of relativistic causal viscous hydrodynamics. 
Eqs. (2. 27-2.29), (2.47-2.49), (2.56), (2.64-2.66), (2.78) and (2.81) represent our main results 
describing the correlation functions for various models of relativistic causal hydrodynamics. 
First we should like to note here that the form of the correlation functions given by 
Eqs.(2. 27-2. 29), (2.47-2.49), (2.56), (2.64-2.66), (2.78) and (2.81) are strikingly similar to 
the correlation functions obtained using relativistic Navier-Stokes theory [3,6]. It should be 
noted that Navier-Stokes equations and the other relativistic causal hydrodynamics theories 
are consistent with the condition that four-divergence of the entropy-current Td^S^ remains 
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Figure 1: (a), (b) and (c) show time evolution of the correlation functions for NS, MIS, DKR 
and JBP hydrodynamics for the initial temperature T^ = 310 MeV . The coefficient of viscosity 
is calculated using 'qoKR = %■ The scaling -q^is = Vjbp = '^I'^^Vokr and t^jvs = '^Vdkr 
ensure that the cross-section remains same in the comparison between the models of hydrodynamics. 
Case (a) corresponds to ^^sf^ = 0.08, while cases (b), (c) depicts for ^^^^ = 0.24, ^^^^ = 1.60 
respectively. In all the figures, solid (red), dotdashed (blue), dotted (green) and dashed (purple) 
curves correspond to the NS, MIS, DKR and JBP Hydrodynamics respectively. 

always positive. The expression for divergence of the entropy four-currents (see Eq.2.17) 
contains tensors like TTny, q^ and scalar IT. Since the derivation of the Onsager coefficients 



7 uses A^^TT 



fiu 



0, Qf^U^ 



and u^TT^'^ 



as constraints on the dissipative fluxes, the 



form of the correlation function (2. 27-2.29), (2.47-2.49), (2.56), (2.64-2.66), (2.78) and (2.81) 
remains similar to the one obtained using the relativistic Navier-Stokes theory. But in case 
of DTT kind of hydrodynamics, it is not clear if divergence of the entropy four-current can 
be expressed directly in terms of scalar product of the viscosity and heat-flux tensors. 

Since forms of the equations of relativistic causal hydrodynamics of MIS, DKR, JBP 
are different, the space-time evolution of the correlation functions can be very different 
from the evolution of the correlations in the Navier-Stokes theory. Next, we have shown 
that when one considers the third order correction [17] to the MIS hydrodynamics, the 
form of correlation function remains to be same as the one obtained using the MIS theory. 
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Figure 2: (a), (b) and (c) show the time evolution of the correlation functions for NS, MIS, DKR 
and JBP hydrodynamics. For Fig. (2a), the ratio of the viscosity to entropy density ^ ~ 0.08 is 
kept same for all the four cases, while figure (2b) and (2c) correspond to cases with ^ = 0.24 and 
^ = 1.60 respectively. Here, solid (red), dotdashed (blue), dotted (green) and dashed(purple) curves 
correspond to the NS, MIS, DKR and JBP Hydrodynamics respectively for all the figures. 



Further, we have considered the case of conformal hydrodynamics when the bulk-viscosity 
is absent. In this case also the form of the Onsager coefficients and the viscous correlations 
remains same as that of MIS theory. Again the temporal evolution of the correlations may 
be different than the MIS theory as the underlying hydrodynamical equations are different. 

As there exists a variety of the relativistic causal hydrodynamics models, it is natural 
to ask the question that which kind of hydrodynamics to be applied for a given physical 
situation. This question is rather hard to answer in general. In order to understand the 
evolution of the correlation functions in some details we have calculated the normalized 
correlation functions given by Eqn.(3.11) for an expanding one-dimensional boost-invariant 
(Bjorken) flow. In figures (1-2), we plot the normalized correlation function C{t) (Eqn.3.11) 
for MIS, JBP, DKR and NS hydrodynamics as a function of time r. The initial values of 
the temperature, Tj = 310 MeV and the viscous stress tt = are kept same for all the 
cases shown in figures (1-2). Each figure has curves representing the four different models 
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of hydrodynamics: the soHd hne describes the correlation function from the Navier-Stokes 
equations, the dotted-dashed curve represents the case of MIS equations, the dotted hne 
describes the DKR case, while the dashed curve describes JBP case. There are two kind of 
comparisons possible between the correlation functions for the different models of relativistic 
viscous hydrodynamics. In one such comparison the energy-independent cross-section [see 
Eqns. (3.16-3.18)] is kept same for all the different versions of the hydrodynamics [19]. 
Following [19,24] one can write the relation between the viscosity coefficient for the different 
models of hydrodynamics as tjokr = %, Vmis = H = Vjbp and rj^s = 0.8436|^. The 
scaling rj^jjg = rjj^p = 9/Wr]£)XR and r^jvs = '^^Vdkr ensure that the cross-section 
remains same in the comparison between the different models of hydrodynamics. In another 
way of comparing the correlation- functions, the ratio of viscosity coefficient to the entropy 
density r]/ s is kept same for the different models of the hydrodynamics, while the transport 
cross-section is varied for the different hydrodynamics models. 

Figure (1) corresponds to the case when the cross-section is kept same for all the 
models of the hydrodynamics. The case(a) corresponds to 2aKfl = 0.08, while r]/ s for 
the other hydrodynamics models is calculated using the scaling. In this case, all the four 
curves overlap with each other. This case corresponds to the nearly ideal hydrodynamics. 
The case (b) and (c) corresponds to ^^ejlr = q 24 and ^^^"^ = 1.6 respectively. One can 
notice here that the correlation function corresponding to the Navier-Stokes hydrodynamics 
consistently gives higher values of the correlation for values of ^cica > 0.08. The DKR curve 
also gives higher value of the correlation in comparison with JBP and MIS hydrodynamics. 
The figure inset in case(b) shows the plots of DKR, JBP and MIS correlation functions in 
T range between 5/m/c to 7/m/c. The reason why the MIS and JBP dynamics follows 
each other is simple as values of parameters /3^ and r^ are same for both the cases and 
dynamics of viscous stress is governed by the Eq.(3.15). In case (c) one can see that 
there is a significant difference between the Navier-Stokes equation and the other causal 
hydrodynamics models. The Navier-Stokes correlation rises for such a high value of the 
viscosity. This is an unphysical behavior. It is well-known that when the Navier-Stokes 
equation is applied to one- dimensional boost-invariant flow, temperature of the flow can 
increase for some time and reach a maximum value before it starts to plummet [26]. 

Figure (2) corresponds to the case when the ratio of the viscosity coefficient to the 
entropy density is kept same for all the four models of hydrodynamics. Case (2a) corresponds 
to the situation when - = 0.08 for all four models of the hydrodynamics. In this case, 
when the viscosity is very low, all the correlations C{t) for all the four models (NS, MIS, 
DKR and JBP) of hydrodynamics overlap. Figure (2b) corresponds to the case when 
^ = 0.24. Here, again the Navier-Stokes hydrodynamics consistently gives higher values of 
the correlation function. The inset diagram in figure (2b) shows the plots of DKR, JBP and 
MIS correlation functions in r range between 5/m/c to 7/m/c. Figure (2c) corresponds to 
case when ^ = 1.6^ here the Navier-Stokes dynamics gives the unphysical behavior. 

The cases studied in figures (1-2) show that Navier-Stokes hydrodynamics gives consis- 
tently higher values for the correlation. The difference between the correlation functions 
obtained using the Navier-Stokes dynamics and of the correlation functions from the other 
causal approaches (MIS, DKR and JBP) increases with increasing values of r]/ s. However, 
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the differences between the correlation-functions of different causal hydrodynamics are not 
very significant. This could be due to the fact that the present example deals with the 
boost-invariant one-dimensional flow. In a non-central heavy-ion collision, the vorticity in 
the flow can play a signiflcant role [27]. As one can notice from Eq.(2.68) that vorticity can 
drive dynamics of the viscous stress in DKR hydrodynamics. However this will require to 
solve hydrodynamical equation in 2 + 1 or 3 + 1 dimensions. This is at present, beyond the 
scope of this work. 

Thus we have shown using one-dimensional flow as an example that evolution of the 
correlation functions of fluctuations can sensitively depend on the model of hydrodynamics. 
The well-separated nature of the correlation-function of the relativistic Navier-Stokes 
equations in comparison with the models of the causal hydrodynamics may allow us, we 
believe, in flguring out which version of the relativistic dissipative hydrodynamics might be 
suitable for describing the heavy-ion collision dynamics. 

5. Conclusions 

We have studied fluctuations in various models of relativistic causal hydrodynamics. We 
have found that the general properties of the dissipative part of the energy-momentum 
tensor due to the viscosity and heat-flux play an important role in determining the Onsager 
coefficients and the correlation functions. We find that all the models of causal relativistic 
hydrodynamics have similar form of the correlation function. In addition, we have applied 
these results to one- dimensional boost-invariant hydrodynamic flow to investigate the 
temporal evolution of the correlation functions for the different hydrodynamics models. We 
believe that our study may help in choosing a proper version of the causal hydrodynamics 
for a given physical situation. 

Note added: After this manuscript was prepared, we have found in Ref [28] on arXiv 
that the authors have applied the fluctuation-dissipation relation to the relativistic viscous 
hydrodynamics with the memory effects. 
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